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In gravitational-wave interferometers, test masses are suspended on thin fibers whicii experience 
considerable tension stress. Sudden microscopic stress release in a suspension fiber, which I call a 
'creep event', would excite motion of the test mass that would be coupled to the interferometer's 
readout. The random test-mass motion due to a time-sequence of creep events is referred to as 'creep 
noise'. In this paper I present an elasto-dynamic calculation for the test-mass motion due to a creep 
event. I show that within a simple suspension model, the main coupling to the optical readout 
occurs via a combination of a "dc" horizontal displacement of the test mass, and excitation of the 
violin and pendulum modes, and not, as was thought previously, via lengthening of the fiber. When 
the creep events occur sufficiently frequently and their statistics is time-independent, the creep noise 
can be well-approximated by a stationary Gaussian random process. I derive the functional form of 
the creep noise spectral density in this limit, with the restrictive assumption that the creep events 
are statistically independent from each other. 



I. INTRODUCTION 

Gravitational-wave interferometers like Laser Interfer- 
ometric Gravitational-wave Observatory (LIGO) in the 
United States of America \\\, VIRGO [5] in Europe, 
and their smaller counterparts GEO600 in Germany [3^ 
and TAMA in Japan [4l, are using super-precise opto- 
mechanical measurements to search for astrophysical 
gravitational waves. After several years of taking sci- 
entific data, LIGO and VIRGO are currently being up- 
graded with improved instrumentation and should again 
become operational in 2015 [5], [6]. LIGO Science Collab- 
oration (LSC) and the VIRGO community are project- 
ing [7] that with the upgraded technology, both inter- 
ferometers will soon be measuring multiple coalescences 
of relativistic compact objects (neutron stars and black 
holes). These projections are based in part on the the- 
oretical predictions for spectral density of the interfer- 
ometers' noise. It is thought that the random processes 
that contribute most of the noise, i.e. the seismic shak- 
ing of the suspensions the thermo- mechanical and 
thermo-refractive fluctuations of the mirror surface [S], 
and the quantum-mechanical fluctuations of the light- 
field coupled to the test-mass motion [10], [11] are well 
understood [T^ . 

One of the dangerous unknowns for the advanced 
gravitational-wave interferometers is a non-Gaussian 
noise from a superposition of transient events in the in- 
strument. In this paper I concentrate on the creep noise, 
which is caused by a superposition of the sudden local- 
ized tension stress releases (creep events) in suspension 
fibers and their end attachments. It has been thought 
that a creep event would couple to the interferometer's 
readout via lengthening of the fiber [TBI. Specifically, 
it was argued that because of the Earth' curvature, the 
laser beam was not strictly perpendicular to the suspen- 
sion fiber, and thus the fiber's lengthening would result 



in some test-mass displacement along the beam. In this 
paper I show that this coupling, while present, is not 
dominant, at least for a simple model where the fiber is 
represented by a cylinder with constant radius. Instead, 
a creep event couples to the interferometer's output pre- 
dominantly through excitation of the pendulum and vi- 
olin modes of the suspension; this coupling is explicitly 
calculated in this work. 

The fact that creep events couple to the transverse vi- 
brational modes of the system is in agreement with the 
experiment of |17j who find a substantial excess noise in 
the transverse motion of a tungsten wire stretched to 20% 
of the break-up stress. Similar excess noise in steel wires 
was observed by [18]. However, the results in [17 were 
not confirmed by |19) who did not observe any excess 
noise in the motion of the stressed tungsten wire. More- 
over, it is far from obvious that the processes responsible 
for the creep events in metallic fibers [2D] will be op- 
erating in the fused silica suspension fibers such as the 
ones that are currently used in GEO600 and that will be 
used in the advanced LIGO, VIRGO, and KAGRA sus- 
pensions [H]. Two experiments with fused silica fibers 
have been performed by |22j and |19j : in both experi- 
ments no excess noise was discovered near violin resonant 
frequency of the fiber. In a more recent work [23], the 
motion of a test mass was monitored in GEO600, where 
the fused silica suspension fibers were used. The motion 
near the violin-mode frequency was entirely consistent 
with that of the thermally excited violin mode. There- 
fore, currently there is no experimental evidence that the 
creep excess noise in the future advanced gravitational- 
wave interferometers will pose a serious problem. How- 
ever, there are at least two reasons to keep investigating 
the creep noise: (1) the measurements in [52], [H], and 
|23) have all been performed at frequencies from several 
hundreds to thousands of Hz, where the noise of ground- 
based interferometers is strongly dominated by the quan- 



2 



turn shot noise, while the danger from creep noise exists 
at much lower frequencies, in the same region of tens of 
Hz where the shot noise is unimportant, and (2) the main 
source of creep noise may well be not inside the fused sil- 
ica suspension fibers, but inside other carrying parts of 
the system like the bond between the test masses and the 
"ears" that are supporting them (Riccardo DeSalvo and 
Norna Robertson, private communications). It is thus 
important to understand how a creep event inside the 
suspension couples to the horizontal motion of the test- 
mass, as well as the frequency dependence of the noise 
generated by a multitude of the creep events. This paper 
lays a theoretical foundation for addressing these issues. 

The plan of the paper is as follows. In section 2, I 
present a convenient reciprocity relationship for linear 
elasto-dynamic systems. In section 3, I use this relation- 
ship to derive the interferometer's response to a creep 
event, as a function of location of the stress release in 
the fiber. In section 4, I derive the functional form of 
the creep noise spectral density, in the limit where the 
creep noise can be treated as a stationary Gaussian ran- 
dom process. A brief discussion of the future work is 
presented in section 5. 



II. ELASTODYNAMICS AND RECIPROCITY 
THEOREM 

During the creep event, the stress changes suddenly in 
some small volume of the fiber. But how does this affect 
the motion of the test mass? At first glance, this seems 
like a formidable problem in elastodynamics. However, 
it turns out that solving the reciprocal problem is suffi- 
cient. Namely, one should in a mental experiment apply a 
sudden force to the test mass and compute the motion of 
the fiber at the location where the creep event originated. 
The solution of the reciprocal problem leads directly to 
the solution of the original problem. Reciprocity rela- 
tions have been thoroughly studied in elastodynamics; 
see e.g. [23] for a comprehensive review. Here I will use 
the following formulation of the reciprocity theorem: 

Consider an elastodynamic system initially at rest, 
that is being driven by a distributed force, with the force- 
per- volume given by 



F{r,t) = f{r)x{t), 

where is some function that is non-zero only 
for t > 0. Consider also a readout variable 



Xit) 



(1) 



(2) 



where ^ (r , t) is the displacement from rest at location r 
and time t. Both the applied forces and displacement 
are assumed to be small, so that a linear approximation 
of elastodynamics holds. The reciprocity theorem states 



that if in the pair of Equations ([T]) and (|2| the form- 
factors f{r) and g{r) are interchanged then the read- 
out variable X{t) remains the same. In other words, the 
input-output dynamical system is invariant with respect 
to / and g interchange, with being the input and 
X{t) being the output. 

The proof of the theorem is as follows. Let ^„(r) be 
the normal modes of the system, with proper angular 
frequencies a;„. A displacement field ^ can then be de- 
composed into a series 



^(r,t) = i;„a„(t)^„(r^. 



(3) 



The mode amplitudes a„ (t) are the new dynamical coor- 
dinates. The Lagrangian of the unforced system is given 
by 



2 21 



(4) 



where m„ is the effective mass of the n'th mode. Ex- 
ternal forcing of Eq. ([T]) is introduced via an additional 
interaction Langangian term 

A„t= J d?rF{r,t) ■ i{r,t), (5) 

which, in terms of the coordinates a„ can be rewritten as 



(6) 



(7) 



The full Langrangian allows us to immediately obtain the 
equations of motion: 



Lint — '^nO-nfnX{t)- 

Here /„ are constants given by 



Therefore, 



a„(i) 



xit) 



In 



-Xn{t), 



(8) 



(9) 



where Xn(t) is the solution to the forced harmonic oscil- 
lator problem 



d'^Xn , 2 
-^+-nXn 



xit) 



(10) 



with the initial condition Xn{0) — Xn(0) = 0. The read- 
out variable in Eq. ([2| can then be written as 



X{t) 



9nfn /,s 
S„ Xn(t), 



where gn is defined similarly to /„: 



(11) 



(12) 



The readout variable X{t) is invariant with respect to 
the interchange of / and g. Q. E. D. 
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III. TEST-MASS RESPONSE TO A SINGLE 
CREEP EVENT 

General considerations 

A creep event happens when a minute section of the 
suspension fiber refuses to support its full share of the 
tension stress. What exactly happens microscopically is 
poorly known, but a simple model will suffice for mod- 
elling of the elastodynamical behavior. Let us assume 
that a small fiber element of volume V suddenly does not 
support any elastic stress Tij. I now consider a slightly- 
reduced elastic system, namely the original one with the 
small volume element V taken out. This slightly-reduced 
system experiences a sudden force applied to the bound- 
ary of the volume element V, so that the boundary sur- 
face element dS, assumed to be directed outside of the 
volume, experiences the force 



dF = -Tij{dS ■ ej)ei, 



(13) 



where are the unit vectors along the coordinate axes, 
and the summation over the dummy indices is assumed. 
I would like to evaluate the test-mass displacement 



X{t) under the action of the force in Eq. (13 1 that is 
switched on at i = (this situations is somewhat similar 
physically to the excitation of magnetar motion as a re- 
sult of sudden reconfiguration of the magnetosphere dur- 
ing a giant magnetar flare; see |25j). By the reciprocity 
theorem from the previous section, this is equivalent to 
acting with the suddenly switched on force on the test 
mass, directed along the laser beam: 



Ft, 



At) - Bit) 



(14) 



where Q{t) is the Heavyside function ^B]. One then has 
to find the response of the slightly- reduced elastic subsys- 
tem to this force, and in particular that of the reciprocal 
readout variable ^readout (t) that is dictated by the func- 
tional form of the force in Eq. ( 13 1 



cadout 



(15) 



where the integration domain is the boundary of the vol- 
ume V . It is obvious that for sufficiently small [27] vol- 
ume V the response of the slightly-reduced system is the 
same as that of the full system, and from hereon I shall 
make no distinction between the two. 

By Gauss' theorem, for small V the above equation 
can be written as 



readout 



-VT,, 



'dxj 



(16) 



where the strain d^i/dxj is evaluated at the location of 
the creep event. To sum up: by finding the response 

of ^readout from Eq. (fl6]) tO the force Ftcstmass = 6(t) 



applied at the test mass along the direction of the laser 
beam, one finds the test-mass displacement in response 
to the creep event. 

It is convenient and instructive to work in the Fourier 
domain 



Ft. 



At) 



dujFt, 



{uj)e 



iu)t 



(17) 



For the force given by Eq. (14), the force Fourier compo- 
nent is given by 

^lcstmass(w) = lim | — — — I , (18) 

e^o+ 2nT[uj — le) J 

where positive e serves to avoid the singularity in 
Eq. ([it]); the hmit e ^- 0- would give Ftcstmass(i) = 



Q{-t). 

The test mass horizontal displacement ^tostmass under 
the action of the applied force is given by 



Xt. 



^Uj) 



M 



Z[{u! — ij{uj)] 



(19) 



where the mechanical impedance Z{itj), as derived in the 
Appendix, is given by [28^ 



Z[lo] 



TT^COt 



— 



(20) 



Here, M is the mass of the test mass, ujp = \/g/l 
is the pendulum angular frequency of the test mass, 
(jjs = "KyJ M / [NrnjiUp is the fundamental violin mode 
angular frequency measured when the test-mass is fixed 
in space, and I, m, and N are the length, the mass, 
and the number of the strings on which the test mass 
is suspended. The small positive 7(w) <C w inserted into 



Eq. ( 19 ) represent damping. Mathematically, this dis- 



places the poles of Z{uj) into the upper half of the com- 
plex w-plane J29]. These poles represent the frequencies 
of normal modes of the suspension; their imaginary parts 
equal the rate of exponential decay of their amplitudes. 
The actual values of 7(0;) are only important near the 
normal-mode frequencies and can be measured experi- 
mentally. 

If no damping is present, the poles of Z[lS\, i.e. the 
normal-mode frequencies, are given approximately by 



~ iwp — ±\/ g/l 
for the pendulum mode, and 



± 



JCJ, -f (-1) 



(21) 



(22) 



for the i — 1,2, ... violin modes. Let us introduce non- 
zero 7p and 7„j which are the damping rates of the pen- 
dulum and violin modes, respectively. The impedance 
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can be expanded as follows: 



Z{u! — ij) 



1 



+ 



(23) 



Substituting Eqs (20) and (18) into Eq. (19), and evalu- 



ating the inverse Fourier transform, I get 



1 



(24) 



[1 - cos(a;i,jt)e-'^"'*] 



As expected, the system's response to a sudden force is a 
constant displacement added to damped oscillations due 
to excited pendulum and violin modes. 

Lets take stock, and sum up what has been done so 
far. In a reciprocal problem, one has to find out the 
induced shear at the location of the creep-event source 
when a sudden force F{t) = Q{t) is applied to the test 
mass. In this subsection I have done part of the problem, 
i.e. I found the test-mass displacement under the action 
of the said force. To make further progress, I need to 
choose a particular model for the suspension fiber itself. 
In the next subsection I consider one of the particular 
cases that can be dealt with analytically. The treatment 
of complicated geometries is left for future work. 



Example: creep event in a cylindrical fiber with 
constant cross-section 

Let us consider in detail the case where the creep events 
occur inside a cyllindrical fiber of constant crossection 
that is rigidly attached at the top to a suspension iso- 
lation plate and at the bottom to the test mass. It is 
assumed here that the allowed test-mass motion is a par- 
allel translation but not rotation, as is the case when 
four suspension fibers are used. The dominant part of 
the stress in the fiber is 



T — 



Mg 



(25) 



where N is the number of suspension fibers and r is the 
radius of the fibers' horizontal crossection. The readout 



variable in Eq. ( 16 ) is given by 



readout 



MgVd^ 



(26) 



where z is the vertical coordinate along the fiber. I choose 
z = at the fiber's top and z = 1 ai the fiber's bottom. 

Let XcUcZ be the spatial coordinates of the source 
of a creep event, where XcUc = corresponds to the 
location of fiber's axis at z, and Xc is measured along the 



laser beam direction. The vertical strain induced by the 
fiber's motion is given by (see, e.g., Chapter 11 of [30] ) 



dz 



(27) 



where 77(2) be horizontal displacement of the fiber. 
Therefore, the readout variable is 



X 



readout 



MgVxc d^T] 



(28) 



Let us now find the fiber motion ri{z,t). Its dynamical 
equation of motion is given by (see, e.g., chapter 12 of 
[Sb] or [21]) 



dz"^ 



A 



dz^ 



(29) 



Here Cs — \/ Mgl/{Nnn) is the velocity of the tension 
wave in a fiber, and A is the characteristic bending length 
given by 



A = 



1 ttEN 



1 



Mg 



^-1/2, 



(30) 



Here r is the radius of the fiber, E is the Young modulus, 
and ^oz.z is the intitial stretch factor of the fiber under the 
loading force Mg/N of the test mass. Advanced LIGO 
will use the fused-silica fibers with the following param- 
eters: r = 2 X lO-^m, M = 40kg, E = 72GPa, I = 0.6m, 
and = 4 fibers per test mass. This parameters produce 
the bending length A ~ 0.001m < /. 

The periodic solutions of the homogenious equation 



( 29 ) can be written as 



'q{z,t) cx e"^*eP^ 



where 



1± Wl 



2A2 



(31) 



(32) 



For frequencies of interest, wA/c^ <C 1, and thus the solu- 
tions feature two physically distinct branches, p = ±1/A 
and p = ibiw/cs. The former branch represents evanes- 
cent quasi-static bending perturbations that will be large 
near the fiber's attachment points, while the latter rep- 
resent tension waves in a fiber. The boundary conditions 
77(0) = 77'(0) = iil) and rjil) = Xtestmass, together with 
A <C ^, determine the full solution for the fiber: 



77(z,t) = Be^'^*[77bc„d(^) + ?7wavo(z)], (33) 



where 



??bcnd(z) = fcA |e ^ — [cos(fcZ) -(- fcAsin(fcZ)]e ^ | 



(34) 
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Coupling of the creep event to test-mass displaceme 




FIG. 1. The quantity C{z,ij) that characterizes the coupUng strength of the creep event to the test-mass horizontal displace- 
ment, is plotted as a function of the creep-event's distance from the top of the fiber, at a frequency of the fundamental violin 
mode. The coupling is strongly peaked within the bending length A from the top. Similar "bending" peak occurs at the bottom 
of the suspension fiber (not shown). The parameters for the plot are those for advanced LIGO suspension fused silica fibers: 
E = 72GPa, r = 2 X W'^m, I = 0.6m, M = 40kg, iV = 4. 



and 

Vwavciz) — sin(fcz) — fcAcos(A;2;). (35) 
Here k — io/cg, and the amplitude B is given by 

B ~ Xtcstmass(w)/ Sm{kl) 



^tcstmass(w)/ Sm{T: UJ / UJ s) . 



(36) 



Therefore, in the Fourier domain the readout variable is 
given by 



^readout f'^) 



lim 



1 



gVxc 



0+ 2T:i{uj — ie) Nnr^ 
1 



TT— COS 


TT — 


— cj^ sin 


TT — 











[a ^77bond(2,a;) - fc^77wavo(2,a;)] , (37) 



where r]\jcnd{z , uj) and ?7wavc(-z, i^) are given by Eqs. (34) 
and ( 35 ) , respectively. In the time-domain, putting in all 



the damping terms, I get 

^readout (0 = Cp{xc, z) [l - cos(a;pt)e"''''*] ^ (38) 
-_xCvj[xc., z) [l - cos(tj,„jt)e^'^''^*] 



V'OO 



where 



Cpi^Xc, z^ 



gVxc 



— ^ s-i 2 2 



j1 P I 



and 



2gVxc{-iy+' 




(39) 



(40) 



In the expressions above I used A ^ I and Wp <C Wg. I 
remind the reader that the coordinates of the 

location of the creep event with the effective volume V , 



and that Eq. (38) gives the test-mass response to such 
an event. The geometric nature of the prefactors in Eqs 
(39) and (40) is apparent once one recalls g/uj^ — I, and 



It is instructive to compute a numerical example. In 
the expression (38) above, consider values Xc = r, V = 
nm"^, and z = (i.e., a formation of a nanometer-size hole 
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at the top edge of the fiber). The displacement that one 
then gets at a pendulum frequency is or order 10~^^m. 

It is worthwhile to have another look at the right-hand 
side of Eq. (37 1. The part of the equation in square 
brackets, 



C(Z,W)=A ^77bond(2:, W) - fc^77« 



(41) 



determines the z-dependence of the coupling of the creep 
event to the horizontal motion of the test mass. The func- 
tion C{z) is plotted in Fig. 1, for uj = uJs (i.e., the funda- 
mental violin mode) . Fiducial parameters that were used 
in making the plot are specified in the figure's header. 
The function peaks very strongly within A from the at- 
tachment ends of the fiber; there C{z) ~ fc/A is domi- 
nated by the ?7bend- Away from the attachment points, 
the coupling is dominated by the jy^avc part of the solu- 
tion and C{z) ~ fc^. 

In this subsection's model the creep events are assumed 
to be triggered homogeneously in the suspension fibers 
|32j . Thus the creep events have only ~ X/l chance to 
be triggered within A from the attachment points. How- 
ever, they have individually much larger impact [by a 
factor of l/(fcA)] on the test-mass motion then those 
ones originating away from the attachments. It fol- 
lows that the creep events originating within the bend- 
ing regions near attachment points contribute most of 
the creep noise; their contribution is greater by a factor 
of ~ k~-^l~^X~^ = {1/'k'^\){uJs/(jj) than that of the creep 
events away from the attachment points. This is studied 
in section 4. 



The case of non-orthogonal laser beam and 
suspension fiber 

Let us now consider the case where a laser beam is 
inclined by a small angle /3 with respect to the horizon- 
tal direction. This is an inevitable effect because of the 
spherical shape of the equipotential surface on which the 
test masses in the same arm are located; for 4km arm 
/? ~ 3 X 10~^ radians. It is this misalignment that was 
previously thought to be the major source of the creep 
noise jl6) ; we treat this mechanism within the formalism 
developed in the previous section. The reciprocal force 
applied at the test mass has now a vertical component 



Fvert(i) = -/3e(i) 

that causes the vertical test-mass motion 



[1 - cos(wvortOe" 



-Tvorttl 



(42) 



(43) 



where the negative sign corresponds to the upward mo- 
tion. Here we take into account only one vertical suspen- 
sion mode with the angular frequency 



Wvert 



Ml 



(44) 



and the damping rate 7vert; the higher-order vertical 
modes are at much higher frequency and have a much 
weaker coupling to the sudden force. The readout vari- 



able from Eq. ( 26 ) is given by 



vcrtrcadout 



£;(iV7rr2)2 
In the Fourier domain, 



(45) 



X 



vcrtrcadout 



1 



Nnr^l 2Triu}[u 



{uj - i7vcrt)^ 



(46) 

It is instructive to compare the amplitude of the vertical 
mode to the amplitude of the pendulum mode excited 
by the creep even near the attachment point, as inferred 



from Eq. ( 38 ) . Their ratio is approximately given by 



vertical 
pendulum 



few X 10 ^. 



(47) 



It is the smallness of this ratio that makes the contri- 
bution to the creep noise from the creep-induced fiber- 
lengthening be subdominant relative to the direct hori- 
zontal coupling, in most of the LIGO band. 



IV. CREEP NOISE IN THE STATIONARY LIMIT 

Consider now a situation where multiple creep events 
are triggered in sequence,. According to the Central 
Limit Theorem, if the events occur sufficiently frequently, 
their superposition produces a random Gaussian noise in 
the test-mass motion. The response of the test-mass to a 
single creep event can be written as X{d, i — ig) where Iq 
is the time when the creep event is triggered, and a is the 
set of parameters characterizing the event (location in the 
fiber, effective volume, etc.). In what follows we assume 
that the creep events are statistically independent from 
each other and that the creep-event parameters sample 
some well-defined probability-distribution function. This 
assumption is known not to hold in some systems that 
exhibit so-called "crackle noise" [33] , and will be relaxed 
in future work. If the probability density distribution 
P(d) is time- invariant, then the creep noise is stationary 
and has a spectral density given by 



Sxif) = 8tt^R J dd\X{d,uj)\^P{d), 



(48) 



where R is the rate of the creep events. The dd implies a 
multi-dimensional integral over the parameter space of a. 
Evaluating this expression for the model of the cylindrical 
fiber, we get 
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frequency (Hz; 



FIG. 2. Two amplitudes of the creep noise are plotted: that due to the direct horizontal coupling (continuous), and that due 
to the vertical lengthening of the suspension fibers (dashed). While the units on the j/-axis are arbitrary, the ratio of the two 
contributions depends on the elastodynamics only and is robust. One can see that the horizontal coupling makes dominant 
contribution everywhere except in a narrow band near the vertical resonance of the last stage of the suspension. The parameters 
for the plot are those for advanced LIGO suspension fused silica fibers: E = 72GPa, r = 2 x lO'^m, I = 0.6m, M = 40kg, 
= 4. For this plot, the Q-factor of all the modes is taken to be 10'^; this choice affects the height of the sharp peaks in the 
figure. Realistic Q values will be several orders of magnitude higher. 



where 







TT COS 









uip — 10^ sin I TT- 



(50) 



and 

G(fc) = -'-{[l + cos2(/cO] + fcA[2sin(2fcO + fc/]}. (51) 
At 

Here (y^) is the ensemble average of the (volume)^ of 
the creep events in the system. Naturally, this quan- 
tity is meaningful only in our simple model for the creep 
events, however a term like this, representing the mean 
of the squared intensity of the creep events, is expected 
in any generic model for the localised creep events. The 
first term in square brakets on the right-hand side of the 
above equations is due to creep events generated near 
the attachment points, while the second term is that due 
to creep events generated in the fiber's bulk; in both of 
these the terms of order fcA have been neglected. It is 



clear that the noise is strongly dominated by the creep 
events near the attachment point. 

The spectral sensity of noise due to the vertical length- 
ening of the fibers is given by 



*S'xvcrt(y) — 

where 



w2 



(53) 



The plots for ^ Sx (/) and \/5'xvGrt(/) are shown in Fig. 2. 
While the scale on the vertical axis of these plots is arbi- 
trary, since R{V'^) is unknown, the spectral density shape 
and the relative contribution of the two noises are fixed. 
We observe that the direct horizontal coupling induces 
greater creep noise than the vertical motion, at all fre- 
quencies except at a narrow band around / — CLivcrt/(27r). 
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V. DISCUSSION 

In this paper I have provided an elastodynamic calcula- 
tion of the interferometer's response to a creep event, and 
found the functional form of the creep noise in the sta- 
tionary limit. A simple model where the fiber was mod- 
eled as a cylinder of constant radius was considered in 
detail, since this allowed me to obtain analytical expres- 
sions for the test-mass response. Two interrelated quali- 
tative features of this model are worth noting: (1) Creep 
events near the fiber's ends receive a much stronger test- 
mass response in the LIGO band than those at the cen- 
ter of the fiber, and contribute the majority of the creep 
noise, and (2) the dominant coupling to the inteferome- 
ter's readout is via excitation of the violin and pendulum 
modes of the suspension, and not via the lenthening of 
the fiber. I should caution though that these conclusions 
may not hold in a fiber with a more complex dependence 
of the cross-sectional radius r on the height z. In partic- 
ular, the fibers in advanced LIGO suspensions are made 
significantly thicker near the end points, in order to min- 
imize the suspension thermal noise. This thickening will 
reduce the local tension stress, thus reducing both the 
coupling of a creep event to the test mass motion and 
the likelyhood of a creep event to occur. 

In a simple model for the creep noise, I have assumed 
that the creep events are triggered homogeneously in the 
suspention fiber. This may not be the case. The creep 
events may be triggered preferentially (1) at the loca- 
tions where the fiber is welded to the test mass or the 
upper suspension plate, although this is not very likely 
since at the weld the fiber is much thicker than in its 
center (r = 1.5 x 10~^m), so the tension is small (Norna 
Robertson, private communications), or (2) near the lo- 
cations where the ears that support the test mass are 
bonded to it. The bonding material is non-metallic and 
non-glassy and is a potential source of problems (Ric- 
cardo DeSalvo, private communications). In future work 
I plan to explore the spatial distribution of the expected 
creep events, as well as relax the assumption of their sta- 
tistical independence. I plan to also deal with the issues 
of non-Gaussianity of the creep-event triggers; it can pre- 
sumably can be mitigated by considering the output of 
several independent interferometers. 

Some comfort for the advanced interferometers can be 
derived from the fact that experiments [22] , [19] , and [23] 
have not observed any influence of the creep noise on the 
violin-mode motion. We note, however, that all of the 
measurements in question have searched for the creep 
noise at high frequencies corresponding to the resonant 
frequencies of violin modes, from several hundred to sev- 
eral thousand Hz. If the creep events are statistically in- 
dependent from each other, then the expected creep noise 
is red, with ^ Sx if) oc except near resonances; see 



the suspension thermal noise in the case where damp- 
ing of the fiber's motion is structural (see, e.g., [34]). 
Therefore, one may argue that since no creep noise that 
exceeds the suspension thermal noise is observed at high 
frequencies, none is expected to exceed the suspension 
thermal noise at low frequencies as well. This argument, 
however, relies on a very simple model for the creep noise 
that was developed in section 4, and in particular it re- 
lies on the creep events being statistically independent. 
This assumption does not hold in many systems that re- 
lease their free energy via spontaneous acoustic emission 
events (known as the "crackle noise"), see [33] and refer- 
ences therein. Thus further experimental and theoretical 
work is warranted for the low-freqyency domain. 

I thank Vladimir Braginsky for impressing on me the 
importance of non-stationary noise in many a conversa- 
tion that we've had from the first time we met in Septem- 
ber of 1994. I thank Misha Gorodetsky for his insightful 
comments given as part of the internal LSC review for 
this article. I thank Riccardo DeSalvo, Eric Gustafson, 
Jan Harms, Norna Robertson, and Peter Saulson for com- 
ments on the earlier version of this manuscript. I thank 
Rana Adhikari, Yanbei Chen, Kip Thorne, and Sergei 
Vyatchanin for useful discussions. Finally, I thank Sarah 
Levin for proofreading parts of the manuscript before 
submission. The calculations in this paper were com- 
pleted during a visit to LIGO laboratory at Caltech. 
The research was supported by the Australian Research 
Council Future Fellowship. 



Eqs. (49) and (|50[). This is the same scaling as that for 
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Appendix: Response of the suspension to a periodic 
force applied at the test mass 

Here we provide a quick derivation; similar derivations 
for more complicated situations when the test-mass tilt 
is allowed is given in e.g. Appendix A of [35]. Suppose a 
periodic force 

F = Foe'"^* (54) 
is acting on the test mass and induces its periodic motion 

^testmass = XqC^'^* . (55) 

The fiber's motion is given by 



_ sin(fcz) t 
'?(,2:,t) — — 7— T-Aoe . 
sm(K«) 



(56) 



Here k = uj/cg = (Tr/l)uj/ujs is the tension- wave vector. 
The horizontal component of back-reaction tension force 
acting on the test-mass is 



^fiber = i^/oe'"* = -Mgkcoi{kl)X^e"^K 
The second Newton's law gives 

Fo + Ffo ^ -Mu^Xq. 



Substituting the Eq. (|57|) above we get 
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which is equivalent to Eq. ( 20 1 in the text 



